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O^'. Abstract 
O 

- 

The innermost stable circular orbits (ISCO) of coalescing neutron star- 
black hole binary are studied taking into account both the tidal and relativistic 
^ ' effects. We adopt the generalized pseudo-Newtonian potential to mimic the 

' general relativistic effects of gravitation. It is found that the separation of 

the neutron star and the black hole at the innermost stable circular orbit 
is greater than that obtained by using either the Newtonian potential (i.e. 
the case in which only the tidal interaction is included.) or the second post- 
Newtonian equation of motion of point mass (i.e. the case in which the effect 
of general relativity is taken into account but the tidal force is neglected.). In 
equal mass cases, it is found that for a /mi ^ 3.5 with a and mi being the 
mean radius and the mass of the neutron star, the tidal effect dominates the 
stability of the binary system while for a/mi ^ 3.5, the relativistic effect, i.e. 
the fact that the interaction potential has unstable orbit, does. The effect of 
the circulation is also studied by comparing the Roche ellipsoids (REs) with 
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the irrotational Roche-Riemann ellipsoids (IRREs). In the IRREs case, which 
are believed to be the case of coalescing binary neutron stars, it is found that 
in the equal mass binary case, the orbital frequency is 495Hz at the ISCO. 

PACS number(s): 04.30.Db 
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I. INTRODUCTION 



The laser interferometers like LIGO 0, VIRGO §, GEO § and TAMA § are currently 
constructed and the detection of gravitational waves is expected at the end of this century. 
One of the most important sources of gravitational waves for these detectors is coalescing 
binary compact stars such as NS-NS, NS-BH and BH-BH binaries. Each mass, each spin 
and the distance of the binary can be determined by applying the matched filter techniques 
^ to the gravitational wave form of the so-called last three minutes of the binary 0. In 
the end of the last three minutes two compact stars coalesce and nonlinear character of the 
gravity and the tidal effects become important, which will be the most exciting part of the 
coalescing event. 

It is known from the study of orbits of a test particle in the Schwarzschild metric that 
the innermost stable circular orbit (ISCO) exists at r = 6M. For binary cases, Kidder, Will, 
& Wiseman |^ investigated a point-mass binary using second post-Newtonian equation of 
motion, and found that the ISCO of the comparable mass binary is at Vgr ~ 7Mtot where 
Tgr is the separation of the binary in the Schwarzschild like radial coordinate. 

As for the tidal effects, Chandrasekhar 0] studied the Roche limit by using Newtonian 
gravity and treating binary systems as incompressible homogeneous ellipsoids. He found that 
the stars are tidally disrupted before contact at rj = (2. 25 Mtot / pY^^ for equal mass binary 
where p is the constant density of ellipsoids. For the typical neutron star with Mtot = 2.8Mq 
and p = 1 X lO^^g/cm'^, becomes 5.6Mtot- 

Recently Lai, Rasio, & Shapiro |^ have investigated the binary consist of finite size 
compressible stars using approximate equilibria. In a series of their papers, they took into 
account the effect of the quadrupole order deviation of stars and found that the hydrody- 
namic instability occurs at r/^ = 6 ~ 1 .5 Mtot for ^ = 0.5 polytropic neutron star with the 
radius Rq = 2. b Mtot- 

In order to know at what radius the final merging phase begins, all three effects, that 
is, the general relativistic, tidal and hydrodynamic effects should be taken into account 
simultaneously since r^^, rt and rh have similar values of ~ 7 Mtot- For this purpose we must 
solve the fully general relativistic equations, which is a difficult 3D problem in numerical 
relativity [|ll|] although the partly general relativistic results have already been presented 
|r^. To know the qualitative feature of the ISCOs we will use various approximations 
to calculate equilibria of the binary in this paper. The reward is that all the calculations 
can be done analytically so that our results will contribute some physical aspects to the final 
understanding of the ISCO. 

In this paper, we will solve the Roche problem |^ as a model of the binary that consists 
of a finite size star and a point-like gravity source. A gravitational potential which has 
unstable circular orbits will be used as an interaction potential of the binary. Specifically, 
we generalize the pseudo-Newtonian potential proposed by Paczyhsky & Wiita |T^ to mimic 
the general relativistic effects. We think that this model mimics a neutron star-black hole 
binary, and describes its qualitative behavior. 

This paper is organized as follows. In §2 the basic equations necessary for constructing 
equilibrium configurations of the Roche ellipsoids (REs) and the Roche-Riemann ellipsoids 
(RREs) |T5| are derived in the case that the interaction potential is general. In §3 circular 
orbits of the neutron star-black hole binary are calculated using the generalized pseudo- 
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Newtonian potential to mimic the effects of general relativity and the results are shown in 
§4. In §5 our results are compared with those using the Newtonian potential and the second 
order post-Newtonian equation of motion . We use the units of c = G = 1 through this 
paper. 



II. GENERALIZED ROCHE-RIEMANN ELLIPSOIDS 

We first regard the black hole as a point particle of mass 1712 and denote the gravitational 
potential by the black hole as V2{r). To mimic the effects of general relativity by V2(r), 
we do not fix the form of V2{r) at the moment. Next we treat the neutron star as an 
incompressible, homogeneous ellipsoid of the semi- axes 01,02 and 03, mass mi and the 
density pi. We treat the self-gravity of the neutron star (Vi) as Newtonian. Although for 
mathematical convenience we assume the incompressibility for the equation of state, the 
effect of the compressibility can be easily taken into account in an approximate way [Q. 
Following Chandrasekhar p[, we call the neutron star the primary and the black hole the 
secondary. 



A. Second Virial Equations 

We use the tensor virial method ^ , and derive the equations necessary for constructing 
equilibrium figures of this system. Choose a coordinate system such that the origin is at the 
center of mass of the primary, the xi-axis points to the center of mass of the secondary, and 
Xs-axis coincides with the direction of the angular velocity of the binary fl. In the frame of 
reference rotating with Q, the Euler equations of the primary are written as 



dui dP d 
dt dxi dxj 



1 I Vmi + vn2 



+ 2pifieii3Mh (2.1) 



where pi,Ui,P and R are the density, the internal velocity, the pressure and the separation 
between the neutron star and the black hole, respectively. 

Let us expand the interaction potential V2 in power series of Xk up to the second order. 
This approximation is justified if R is much larger than oi, 02 and 03. We assume that the 
potential V2 is spherically symmetric so that it depends only on the distance r from the 
center of mass of the secondary as 

V2 = V2{r), (2.2) 

where r is given by 

r = {{R - xi)^ + xl + xiy^\ (2.3) 
The expansion of V2(r) becomes 
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where the subscript denotes the derivatives at the origin of the coordinates. In the case 
of the circular orbit, we have from the force balance at the center 

'"^''^--(^1(1 + ^). (2.5) 



nil + "^2 \ / 

where 6 is the quadrupole term of the interaction potential 
Substituting Eqs. (|2.4|) and ( p.5|) into Eq. (|2.1|) , we have 



dui dP d 
dt dxi dxi 



1 dW, 



I xi 





+ 2piQeii3Ui. (2.6) 



2R \ dr 

Multiplying Xj to Eq.( p.6| ) and integrating over the volume of the primary, we have 
I / piu^x.d'^ = 2T,, + W,, + |^]^ + (1^) j Suh, 



2fieii3 J PiUiXjd^x + 5ijU, (2.7) 



where 



and 



Tij = 2 J Pi'^i'^jd^^ • Kinetic Energy Tensor, (2.8) 
Wij = J pi-^—Xjd'^x. : Potential Energy Tensor, (2.9) 

lij = J piXiXjd^^ : Moment of Inertia Tensor, (2-10) 

11= ( Pd^^. (2.11) 



In Eq.( |2.7] ) there is no terms related to 5. Since it is possible to take the coordinate system 
comoving with the center of mass of the binary system, the term proportional to 5 in Eq. ( p.6|) 
vanishes when we integrate over the volume of the primary. Eq.( |2.7| ) is the basic equation 
to construct the equilibrium figures of the Roche ellipsoids (REs) and the Roche-Riemann 
ellipsoids (RREs) for the general potential V2(r). 



B. Equilibrium Roche-Riemann Sequence 

In this subsection, we show how to construct the equilibrium figures of the RREs. The 
Roche-Riemann ellipsoid is the equilibrium in which the shape of the primary does not 
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change in the rotating frame although the uniform vorticity exists inside the primary. We 
restrict the problem to the simplest case where the uniform vorticity of the primary is parallel 
to the rotation axis, i.e. the primary is the Riemann S-type ellipsoid 

We set the coordinate axes to coincide with the principal axes of the primary. For the 
uniform vorticity (, the internal velocity in the rotating frame is given by 

ui = Q1X2, (2.12) 
U2 = Q2X1, (2-13) 

and 



Us = 0, 



(2.14) 



where 



Qi 



—C 



of + al 



(2.15) 



and 



Q 



2 ~ — 2' 

af + 02 



(2.16) 



For the stationary equilibrium, Eq.( p.7| ) is rewritten as 



QikQjihi + Wii + <^^ + 



2iJ-2j 



1 (dVo 



R \ dr 



(2.17) 



where Qij is not zero only for 



Q12 — Qi) 
Q21 = Q2- 



(2.18) 
(2.19) 



Eq.( |2.17D has only diagonal components as 



QII22 + Wn + \ 



'd^V2 
Qj.2 



Qlh. + W22 + 1^' + ^ 



In + 2nQ2iu = -n. 



/22 - 21]Qi/22 = -n. 



(2.20) 
(2.21) 



and 



(2.22) 
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We assume for simplicity that the gravitational potential of the primary is Newtonian. In 
this case, the potential energy tensor and the moment of inertia tensor of the incompressible, 
homogeneous ellipsoids are calculated as 



Wij = -2TTpiAiI, 



and 



where 



f 

^0 



°° du 



A [aj + u) 



and 



Eliminating H from Eqs.(|]2g)-(|2]2|), we have 



1 + 2 



a\ + 02' 



1 (dVo 



0-1 — ^ 



a 



.2 „2 



^ R\dr 

271 pi (al - a^] Bi3 



and 



1 + 2^^//? + 
af + 02 



0x02 ^ 



1 /Sr, 



R \ dr 



1 



^2 - 



2 „2 



R\dr 

27rpi (03 - ag) B 



23) 



where 

and the following relations are used; 

ofAi - a]Aj = (a,^ - aj) B, 

where 



Bij = aia2a3 



Now from Eq.(2.5) Q is given by 

^2 l+pfdV2^ 



R \ dr 



p 



nil 
m2 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



(2.32) 



Dividing Eq.( p.27D by Eq.( p.28| ), we have the equation to determine the Roche-Riemann 
sequences as 



(l+p)(l + 5)< 


[1 + 2^ 










(l+p)(l + 5)< 


^1 +2a^+a^/K 1 {aj+aif^) J 







^1 ~ '^3) -^13 
i2 — 03) _B23 



. (2.33) 



Using Eqs. (|2.28| ) and (|2.32|) , we can determine the orbital angular velocity i7 by 
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fi2 _ 2(1 + 5) (a2-a2)B 

[(1 + + ^) {1 + 2^/h + [-^Jr 
Note that //j is related to the circulation C as 

C = /u_.a.-c^l = vra,a,(2 + /^)l] 

where 

{Uinertial)l = {Ql — ^)X2 
i^'^inertial) 2 (^ 



^inertial 



{Uinertial)2 — {Q2 + ^)Xi — 



(2.34) 



(2.35) 



(2.36) 



If there is no viscosity inside the primary, the circulation should be conserved from Kelvin's 
circulation theorem. 



C. Total Angular Momentum 

The total energy and the total angular momentum of the binary are the decreasing 
functions of time since the gravitational waves are emitted. If the total angular momentum 
has its minimum at some separation of the binary, we regard this point as the ISCO [|. The 
total angular momentum of our system which is the sum of the orbital and the spin angular 
momentum is given by 



Jtot = mirl^n + m2{R- rem 

1 



5 af + a2 



mi + m2 [ 5 \ + 02 y 



(2.37) 



where 



^Lai, Rasio, Sz Shapiro show in appendix D of that the true minimum point of the total energy 
coincides with that of the total angular momentum. Strictly speaking, if the rotation includes only 
to the quadrupole order, this coincidence fails. However the difference is as small as the numerical 
accuracy ||9[. 
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m2-R 
mi + 1712 



(2.38) 



The first term in the braces of the right hand side of Eq.( 2.37| ) comes from the orbital angular 
momentum of the binary system and the second does from the spin angular momentum of 
the primary. 



III. GENERALIZED PSEUDO-NEWTONIAN POTENTIAL 

There are variety of choices of V2(r) to mimic the general relativistic effects of the grav- 
itation. We generalize the so-called pseudo-Newtonian potential proposed by Paczyhsky & 
Wiita originally. This potential fits the effective potential of the Schwarzschild black 



hole quite well as we will show later. We will use the generalized pseudo- Newtonian potential 
defined by 

V2ir) = (3.1) 



r — 

rpseudo = rs{l + g{p)} , (3.2) 
^ 7.49p _ lO.y 29.3p3 
^^^^ 6(l+p)2 3(l+p)4 + 6(l+p)6' ^"^-"^^ 



2GM 



tot 



C2 



(3.4) 



Mtot = mi+m2, (3.5) 

where p = mi/m2 and g{j)) is the special term to fit the ISCOs of the hybrid second post- 
Newtonian calculations by Kidder, Will, & Wiseman 0. For p = 0, the generalized pseudo- 
Newtonian potential agrees with the pseudo- Newtonian potential proposed by Paczyhsky & 
Wiita ig. 

Fig. 1(a) shows effective potentials (solid lines) and locations of circular orbits (dots) in 
our generalized pseudo-Newtonian potential ( p = & r pseudo = rg) and in the Schwarzschild 
metric. Although by this choice of the parameter {rpseudo = rg), the locations of the ISCOs 
in the generalized pseudo-Newtonian potential agree with those in the Schwarzschild metric, 
the angular momenta at the ISCO are different, that is, the angular momentum in the 
generalized pseudo-Newtonian potential {J pseudo) for p = is (9/8) times larger than that 
in the Schwarzschild metric (Jsch) at the ISCO. Therefore in Fig. 1(a) and (b) we compare 
circular orbits with different angular momentum related as 

Jpseudo = ( o ) -Jsch- (3.6) 



From Fig. 1(b) we see that the radii of the circular orbits of the generalized pseudo-Newtonian 
potential agrees with those of the effective potential around Schwarzschild black hole within 
10% accuracy near the ISCO. This is the reason why we believe that our generahzed pseudo- 
Newtonian potential expresses the effect of general relativity within 10% or so. 
Using Eq.( |3.1| ) and Eq.( |2.32| ), we can rewrite Eq.( p.34| ) as 
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^2 ~ (^1)^23 



12(1 + p)m{R - Tpseudo) 



- 1 



0, (3.7) 



where a is the mean radius of the primary. In the generahzed pseudo-Newtonian case, the 
quadrupole term 6 is written as 



{3R- 

^pseudo) {R ^pseudo) 



I- — ir)2 yJJ^i^ ' pseudo J y-'^ ' pseudo) / 2 

^ = 10 Y""' 3^^ + 

We also have the separation of the binary as 

E 



1 



{R- 

^pseudo J 



(3.8) 



R 



E-2 



pseudo 1 



(3.9) 



where 



E 



:i+p)(i+5) 1+2 



fR + 



+ 



;i + p)(i + 5) 1 + 2 



al + al 



fR + 




al)B 
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al — al)B' 



. (3.10) 
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For the given mass ratio p, mean radius /bara/nii, circulation parameter fji and axial 
ratio a^/ai, we can determine the axial ratio 02/01 by solving Eq. (|3.7|) with Eqs.( |3.8| ) and 
5.91). Using the axial ratios (02/01, 03/01), we are able to calculate the orbital angular 



velocity by Eq.( 2.34|) and the separation of the binary by Eq. (|3.9|) . The total angular 
momentum is calculated by Eq. ( |2.37| ) . Finding the minimum of the total angular momentum, 
we can determine the location of the ISCO. 

When the viscosity inside the primary is so effective that no internal motion exists, = 
in the above equation, i.e. the Roche ellipsoids (REs). While if the primary is in- viscid and 
C = 0, we have fji = —2, i.e. the irrotational Roche- Riemann ellipsoids (IRREs). Note that 
if we substitute the Newtonian potential as an interaction potential, Eqs.( |2.3^^ ) and ( f^.34[ ) 
agree with the equations derived by Chandrasekhar p in the REs (/r = 0) case and those 
by Aizenman |TH] in the RREs case. 



IV. RESULTS 



Kochanek [T^ and Bildsten & Cutler showed that the internal structure of a coalesc- 
ing binary neutron star is the irrotational Roche- Riemann ellipsoids (IRREs). However we 
calculate both the REs and the IRREs for comparison. We show the results in four cases; 
1) the REs with p = 1 or 0.1, 2) the IRREs with p = 1 or 0.1. Results are given in Table. I 
to Table. VI, where Q = Q/y/Wpi represents the normalized orbital angular velocity, and J 
denotes the normalized total angular momentum defined by 

j= "l^ f4 1) 

mim2(r,/Mto,)i/2' ^ ' ' 
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a is the mean radius of the primary defined by 



( 



) 



1/3 



a = 



(4.2) 



In Tables I to VI, f means the point of the ISCO defined in this paper, and if does the 
point of the Roche limit where the Roche limit is defined by the distance of closest approach 
for equilibrium to be possible |^. The values in the parentheses show the power of 10. 



Fig. 2(a) shows J as a function of the normalized separation R/rg. In Fig. 2(a), thin solid, 
dotted and dashed fines are the REs with a/mi being 3, 5 and 8, respectively, while thick 
solid, dotted and dashed lines are the IRREs with a/rrii being 3, 5 and 8, respectively. We 
defined in §§II.C that the location of the ISCO is the minimum point of J. From Fig. 2(a), 
we see that the separations of the binary at the ISCO in the IRREs case are almost the 
same as those in the REs case. 

Fig. 2(b) shows the axial ratios 02/01 and 03/01 as a function of R/rg. The conventions 
of lines are the same as those in Fig.2(a). The relation among the length of the axes is 
Oi > 02 > 03 in the REs, while Oi > 03 > 02 in the IRREs. In the REs, the tidal force 
makes the Oi axis long while it does 02 and 03 axes short. The centrifugal force makes oi 
and 02 axes long. As a result we have Oi > 02 > 03. In the IRREs, in addition to the above 
effects, the Coriolis force caused by the internal motion of the primary exists. This makes 
oi and 02 axes short, which yields oi > 03 > 02 . 

From Fig.2(b), it is found that in the IRREs the star with o/mi = 3 reaches the ISCO at 
the point of 02/01 ~ 0.912 and 03/01 ^ 0.918. On the other hand, the star with a/nii = 8 
terminates when 02/01 ~ 0.760 and 03/01 ~ 0.792. This means that the primary with the 
smaller mean radius reaches the ISCO before the shape of the primary deviates from the 
sphere considerably. This tendency is the same for the REs. 



Fig. 3(a) shows J as a function of R/vg. The conventions are the same as those in Fig.2(a) 
except a/rrii being 3, 5 and 8. From Fig. 3(a), we see that all lines with the mean radii in 
the range of 3 < o/mi < 8 take their minimum at the points near R/vg ~ 3.25. This value 
is almost the same as the ISCO by Kidder, Will, & Wiseman []^ for p = 0.1. This means 
that when p is much less than 1 and the mean radius of the primary is less than 8mi which 
corresponds to 17km for mi = 1.4Mq, the size of the primary has little effect on the ISCO. 

Fig. 3(b) shows the axial ratios O2/01 and 03/01 as a function of R/vg. The conventions 
are the same as those in Fig. 3 (a). We see that the stars with smaller mean radii reach the 
ISCO even when the deviation from the spherical symmetry is very small. Since the binary 
system enters an unstable circular orbit before the primary is tidally deformed, the tidal 
effects are not important. We see that even o/mi ~ 8 the minimum values O2/01 and O3/01 
are not small (~ 0.85). 



A. p = 1 Case 



B. p = 0.1 Case 
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V. DISCUSSIONS 



In this section, we discuss the differences between the case of the generahzed pseudo- 
Newtonian potential and that of the Newtonian potential, and between the REs and the 
IRREs. We will also compare our results with other papers. 

A. The Case of p = 1 

In Fig. 4, the separation Risco/fs is shown as a function of a/nii. Thick lines and 
thin lines represent the case of the generalized pseudo-Newtonian potential and that of 
the Newtonian potential, respectively. Solid lines and dotted lines express the case of the 
IRREs and the REs, respectively. We see that for the Newtonian potential, Risco/^s 
increases in proportion to d/rrii regardless of types of ellipsoids, while for the generalized 
pseudo-Newtonian potential, the behavior of Risco/^s changes around a/rui ~ 3.5. 

For d/rrii ^ 3.5, Riscol'^'s increases in proportion to d/rrii as in the case of the Newtonian 
potential. In this region, the tidal effect dominates the system and the effects of the general 
relativity become less important. In conclusion, for d/mi ^ 3.5, the tidal effect dominates 
the stability of the binary system and for d/mi ^3.5, the binary system is dominated by 
the relativistic effect, i.e. the fact that the interaction potential has an unstable orbit. From 
Fig. 4 it is also found that the location of the ISCOs in the IRREs case is not so different 
from that in the REs case for the same d/mi. 

The orbital frequency of the IRREs at the ISCO for mi = IAMq and d/mi = 5 is 
estimated from Q in Table. IV as 

nisco = 495 [Hz]. (5.1) 
This value is smaller than that in the Newtonian potential case (f^isco°" ~ ^^9 [Hz]). 

B. The Case ot p = 0.1 

If we consider the primary as the neutron star of mass 1.4M0, then, the secondary is 
regarded as the black hole of mass 14M0. Risco/^s is shown as a function of d/rrii in 
Fig. 5. For the Newtonian potential, Risco/^s increases in proportion to d/mi like the case 
of p = 1. For the generalized pseudo-Newtonian potential, in the range of d/mi of Fig.5, 
Risco/^s converges to the value 3.25 obtained by Kidder, Will, & Wiseman |Q. This is 
because the existence of the unstable orbit in the generalized pseudo-Newtonian potential 
influences Risco/fs, and this effect dominates when the radius of the primary is small. 
This is clearer for the small mass ratio p. Therefore if p is small, for the range of the radius 
(3 < d/mi < 8) relevant to the neutron star, the effect of the neutron star's size is very 
small. This comes essentially from the fact that the Newtonian estimate of the Roche radius 
is smaller than the radius of the ISCO. 

One can estimate the orbital frequency of the IRREs at the ISCO for mi = 1.4Mq and 
d/mi = 5 from Q in Table. V as 

^^isco = 187 [Hz]. (5.2) 
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C. Comparison with Other Works 



Kidder, Will, & Wiseman ^ studied the motion of the point-particle binary systems 
using hybrid Schwarzs child second post-Newtonian equation of motion, and obtained the 
separation at the ISCO (rfsco) expressed by 

~ 6 + 7A9ri - 20.8ri^ + 29.3ri^ (5.3) 

Mtot 

where 

V = (5.4) 



Table. VII shows the comparison of Eg. ( |5.3| ) with our results for a /mi = 5. It is found 



that for p = 0.1, the finite size effect is not important because the primary is much lighter 
than the secondary, so that the results are almost the same as [0. On the other hand, when 
p = 1 the finite size of the primary increases the separation at the ISCO by the general 
relativistic and tidal effects. 



Lai, Rasio, & Shapiro [|^, discussed the relativistic effects on the binary system for 
compressible Darwin ellipsoids using a simple approximate model The Darwin ellipsoids 
are the equilibrium configurations constructed by two identical synchronized finite-size stars 
including the mutual tidal interactions. In their approach the effects of general relativity 
and the Newtonian tidal interactions for finite-size compressible stars are combined by hand. 
While we formulated the problem using arbitrary interaction potentials of the secondary for 
the incompressible primary. We adopted the semi-relativistic potential called the generalized 
pseudo-Newtonian potential to mimic the general relativistic effects of gravitation. We 
solved the equilibria of the REs and the IRREs in this potential. Their results and ours 
are compared in Table. VIII, where expresses the minimum separation obtained by Lai, 
Rasio, & Shapiro |jT8|. From this table, we see that both results agree rather well inspite of 



different approaches and approximations. 
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TABLE CAPTIONS 



Table.I. Equilibrium Sequences of the Roche eUipsoids (REs) and the irrotational Roche 
Riemann eUipsoids (IRREs) with p — 1 and a/rrii — 3. 

Table.II. Equihbrium Sequences of the REs and the IRREs with p = 1 and a/mi = 5. 

Table.III. Equihbrium Sequences of the REs and the IRREs with p—1 and a/rrii — S. 

Table.IV. Equilibrium Sequences of the REs and the IRREs with p — 0.1 and a/rrii — 3. 

Table. V. Equilibrium Sequences of the REs and the IRREs with p = 0.1 and a/mi = 5. 

Table. VI. Equilibrium Sequences of the REs and the IRREs with p — 0.1 and a/rrii — 8. 

Tabic. Vll. Comparison of our results (a/mi = 5) with that of Kidder, Will, & Wiseman[7] 
in the cases of p = 1 and 0.1. 

Table. Vlll. Comparison of our results (p = 1) with that of Lai, Rasio, & Shapiro[18] in 
the cases ol a/rrii — 5 and 8. 
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TABLES 



p = 1 



a/mi = 3 



Roche Sequences 


Irrotational Roche-Riemann Sequences 


0,3/01 


09/01 




J 


R/r-s 




('2 /(I I 




> 


.7 




0.950 


0.966 


1.87(-2) 


2.97 


4.75 


0.950 


0.948 


2.69( 


-2) 


2.87 


4.31 












0.918 


0.912 


4.10( 


-2) 


2.86 


3.87 t 


0.912 


0.939 


3.19(-2) 


2.94 


4.13 t 














0.900 


0.929 


3.59(-2) 


2.94 


4.00 


0.900 


0.891 


4.81( 


-2) 


2.86 


3.72 


0.850 


0.889 


5.19(-2) 


2.97 


3.65 


0.850 


0.833 


6.53( 


-2) 


2.88 


3.46 


0.800 


0.847 


6.67(-2) 


3.01 


3.44 


0.800 


0.773 


7.92( 


-2) 


2.92 


3.31 


0.750 


0.801 


8.03(-2) 


3.05 


3.30 


0.750 


0.715 


9.02( 


-2) 


2.95 


3.22 


0.700 


0.754 


9.23(-2) 


3.10 


3.21 


0.700 


0.657 


9.87( 


-2) 


3.00 


3.17 


0.650 


0.704 


1.03(-1) 


3.15 


3.14 


0.650 


0.601 


1.05( 


-1) 


3.04 


3.14 


0.600 


0.652 


l.ll(-l) 


3.20 


3.10 


0.600 


0.547 


1.09( 


-1) 


3.09 


3.12 












0.577 


0.524 


1.11( 


-1) 


3.12 


3.12 X 


0.550 


0.599 


1.17(-1) 


3.26 


3.08 


0.550 


0.496 


1.12( 


-1) 


3.15 


3.12 


0.514 


0.560 


1.20(-1) 


3.30 


3.07 X 














0.500 


0.544 


1.21(-1) 


3.32 


3.07 


0.500 


0.447 


1.12( 


-1) 


3.21 


3.14 


0.450 


0.489 


1.22(-1) 


3.38 


3.09 


0.450 


0.400 


1.11( 


-1) 


3.28 


3.17 


0.400 


0.433 


1.19(-1) 


3.45 


3.12 


0.400 


0.354 


1.07( 


-1) 


3.36 


3.21 


0.350 


0.376 


1.14(-1) 


3.53 


3.18 


0.350 


0.310 


1.02( 


-1) 


3.45 


3.28 


0.300 


0.320 


1.05(-1) 


3.63 


3.27 


0.300 


0.267 


9.33( 


-2) 


3.56 


3.37 


0.250 


0.264 


9.17(-2) 


3.75 


3.41 


0.250 


0.225 


8.23( 


-2) 


3.69 


3.51 


0.200 


0.209 


7.53(-2) 


3.90 


3.62 


0.200 


0.182 


6.83( 


-2) 


3.86 


3.72 



TABLE I. 
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p = 1 



a/mi = 5 



Roche Sequences 


Irrotational Roche-Riemann Sequences 




02/01 




J 


R/rs 


03/01 


02/01 


02 


J 


R/rs 


0.950 


0.967 


1.97( 


-2) 


3.30 


7.24 


0.950 


0.947 


2.90( 


-2) 


3.12 


6.47 


0.900 


0.931 


3.81( 


-2) 


3.18 


5.99 


0.900 


0.890 


5.25( 


-2) 


3.02 


5.48 


0.850 


0.892 


5.56( 


-2) 


3.15 


5.40 


0.850 


0.831 


7.15( 


-2) 


2.99 


5.05 














0.832 


0.810 


7.73( 


-2) 


2.99 


4.95 t 


0.828 


0.874 


6.29( 


-2) 


3.15 


5.22 t 














0.800 


0.850 


7.19( 


-2) 


3.15 


5.04 


0.800 


0.771 


8.67( 


-2) 


2.99 


4.81 


0.750 


0.806 


8.68( 


-2) 


3.17 


4.81 


0.750 


0.711 


9.85( 


-2) 


3.01 


4.66 


0.700 


0.759 


1.00( 


-1) 


3.20 


4.64 


0.700 


0.653 


1.08( 


-1) 


3.04 


4.57 


0.650 


0.709 


1.12( 


-1) 


3.24 


4.53 


0.650 


0.597 


1.14( 


-1) 


3.08 


4.52 


0.600 


0.657 


1.21( 


-1) 


3.28 


4.46 


0.600 


0.543 


1.19( 


-1) 


3.14 


4.50 














0.582 


0.524 


1.20( 


-1) 


3.16 


4.50 X 


0.550 


0.604 


1.28( 


-1) 


3.34 


4.42 


0.550 


0.492 


1.21( 


-1) 


3.20 


4.51 


0.513 


0.564 


1.31( 


-1) 


3.38 


4.41 X 














0.500 


0.549 


1.32( 


-1) 


3.40 


4.41 


0.500 


0.443 


1.21( 


-1) 


3.28 


4.53 


0.450 


0.493 


1.33( 


-1) 


3.48 


4.44 


0.450 


0.396 


1.20( 


-1) 


3.36 


4.58 


0.400 


0.436 


1.30( 


-1) 


3.57 


4.50 


0.400 


0.351 


1.15( 


-1) 


3.47 


4.66 


0.350 


0.378 


1.23( 


-1) 


3.68 


4.60 


0.350 


0.308 


1.09( 


-1) 


3.59 


4.77 


0.300 


0.321 


1.13( 


-1) 


3.81 


4.75 


0.300 


0.265 


9.96( 


-2) 


3.75 


4.94 


0.250 


0.265 


9.85( 


-2) 


3.98 


4.99 


0.250 


0.223 


8.75( 


-2) 


3.94 


5.17 


0.200 


0.209 


8.03( 


-2) 


4.22 


5.34 


0.200 


0.182 


7.23( 


-2) 


4.19 


5.51 
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p = 1 



a/mi = 8 



Roche Sequences 


Irrotational Roche-Riemann Sequences 




02/01 




J 


R/rs 


03/01 


02/01 


02 


J 


R/rs 


0.950 


0.968 


2.02( 


-2) 


3.82 


11.00 


0.950 


0.947 


3.04( 


-2) 


3.56 


9.71 


0.900 


0.932 


3.95( 


-2) 


3.61 


8.98 


0.900 


0.890 


5.53( 


-2) 


3.37 


8.13 


0.850 


0.894 


5.78( 


-2) 


3.54 


8.03 


0.850 


0.830 


7.54( 


-2) 


3.30 


7.44 


0.800 


0.853 


7.50( 


-2) 


3.51 


7.46 


0.800 


0.769 


9.13( 


-2) 


3.28 


7.06 














0.792 


0.760 


9.35( 


-2) 


3.28 


7.02 t 


0.774 


0.830 


8.34( 


-2) 


3.51 


7.24 t 














0.750 


0.809 


9.09( 


-2) 


3.51 


7.07 


0.750 


0.709 


1.04( 


-1) 


3.29 


6.84 


0.700 


0.762 


1.05( 


-1) 


3.53 


6.81 


0.700 


0.651 


1.13( 


-1) 


3.31 


6.70 


0.650 


0.712 


1.17( 


-1) 


3.56 


6.63 


0.650 


0.595 


1.20( 


-1) 


3.36 


6.62 


0.600 


0.661 


1.27( 


-1) 


3.60 


6.51 


0.600 


0.541 


1.24( 


-1) 


3.42 


6.59 














0.584 


0.524 


1.25( 


-1) 


3.44 


6.59 X 


0.550 


0.607 


1.35( 


-1) 


3.66 


6.45 


0.550 


0.490 


1.27( 


-1) 


3.49 


6.60 


0.513 


0.566 


1.38( 


-1) 


3.72 


6.43 X 














0.500 


0.552 


1.39( 


-1) 


3.74 


6.44 


0.500 


0.441 


1.27( 


-1) 


3.59 


6.64 


0.450 


0.495 


1.40( 


-1) 


3.83 


6.48 


0.450 


0.394 


1.25( 


-1) 


3.70 


6.72 


0.400 


0.438 


1.37( 


-1) 


3.94 


6.57 


0.400 


0.349 


1.20( 


-1) 


3.83 


6.85 


0.350 


0.380 


1.30( 


-1) 


4.08 


6.74 


0.350 


0.306 


1.13( 


-1) 


3.99 


7.03 


0.300 


0.322 


1.18( 


-1) 


4.26 


6.99 


0.300 


0.264 


1.04( 


-1) 


4.19 


7.29 


0.250 


0.265 


1.03( 


-1) 


4.49 


7.37 


0.250 


0.223 


9.06( 


-2) 


4.45 


7.66 


0.200 


0.210 


8.33( 


-2) 


4.80 


7.94 


0.200 


0.181 


7.47( 


-2) 


4.78 


8.22 
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p = 0.1 



a/mi = 3 



Roche Sequences 


Irrotational Roche-Riemann Sequences 


as/ai 


02/01 




J 


R/rs 


03/01 


02/01 


02 


J 


R/rs 












0.992 


0.992 


2.42(-3) 


2.70 


3.26 t 


0.990 


0.992 


2.39(-3) 


2.71 


3.27 t 












0.950 


0.959 


9.94(-3) 


2.85 


2.35 


0.950 


0.949 


1.17(-2) 


2.89 


2.27 


0.900 


0.914 


1.80(-2) 


3.02 


2.08 


0.900 


0.896 


2.04(-2) 


3.06 


2.03 


0.850 


0.868 


2.51(-2) 


3.16 


1.96 


0.850 


0.843 


2. 76 (-2) 


3.19 


1.93 


0.800 


0.821 


3.12(-2) 


3.26 


1.88 


0.800 


0.789 


3.35(-2) 


3.29 


1.86 


0.750 


0.772 


3.65(-2) 


3.35 


1.84 


0.750 


0.735 


3.84(-2) 


3.37 


1.82 


0.700 


0.722 


4.10(-2) 


3.42 


1.80 


0.700 


0.682 


4.23(-2) 


3.43 


1.79 


0.650 


0.672 


4.47(-2) 


3.48 


1.78 


0.650 


0.629 


4.53(-2) 


3.48 


1.77 


0.600 


0.621 


4.75(-2) 


3.52 


1.76 


0.600 


0.577 


4.74(-2) 


3.51 


1.76 


0.550 


0.569 


4.93(-2) 


3.56 


1.75 


0.550 


0.526 


4.86(-2) 


3.54 


1.76 












0.525 


0.501 


4.89(-2) 


3.55 


1.76 X 


0.500 


0.517 


5.02(-2) 


3.58 


1.75 t 


0.500 


0.476 


4.90(-2) 


3.55 


1.76 


0.450 


0.465 


5.01(-2) 


3.59 


1.75 


0.450 


0.427 


4.85(-2) 


3.56 


1.76 


0.400 


0.412 


4.89(-2) 


3.58 


1.76 


0.400 


0.379 


4.70(-2) 


3.55 


1.78 


0.350 


0.360 


4.66(-2) 


3.57 


1.78 


0.350 


0.331 


4.46(-2) 


3.54 


1.79 


0.300 


0.308 


4.29(-2) 


3.54 


1.81 


0.300 


0.284 


4.11(-2) 


3.51 


1.82 


0.250 


0.255 


3.80(-2) 


3.49 


1.85 


0.250 


0.238 


3.64(-2) 


3.46 


1.86 


0.200 


0.204 


3.17(-2) 


3.42 


1.92 


0.200 


0.192 


3.05(-2) 


3.40 


1.93 



TABLE IV. 
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p = 0.1 



a/mi = 5 



Roche Sequences 


Irrotational Roche-Riemann Sequences 




02/01 




J 


R/rs 


03/01 


02/01 


02 


J 


R/rs 














0.962 


0.961 


1.10(-2) 


2.71 


3.28 t 


0.953 


0.962 


1.07( 


-2) 


2.71 


3.30 t 












0.950 


0.960 


1.14( 


-2) 


2.71 


3.25 


0.950 


0.949 


1.38(-2) 


2.71 


3.09 


0.900 


0.917 


2.12( 


-2) 


2.75 


2.79 


0.900 


0.896 


2.47(-2) 


2.75 


2.70 


0.850 


0.872 


3.00( 


-2) 


2.79 


2.58 


0.850 


0.841 


3.38(-2) 


2.80 


2.51 


0.800 


0.825 


3.78( 


-2) 


2.83 


2.45 


0.800 


0.787 


4.14(-2) 


2.84 


2.41 


0.750 


0.777 


4.47( 


-2) 


2.87 


2.37 


0.750 


0.732 


4.76(-2) 


2.88 


2.34 


0.700 


0.728 


5.06( 


-2) 


2.91 


2.31 


0.700 


0.678 


5.17(-2) 


2.90 


2.30 


0.650 


0.677 


5.55( 


-2) 


2.94 


2.27 


0.650 


0.624 


5.64(-2) 


2.93 


2.26 


0.600 


0.626 


5.92( 


-2) 


2.97 


2.24 


0.600 


0.572 


5.90(-2) 


2.95 


2.24 


0.550 


0.574 


6.18( 


-2) 


2.99 


2.23 


0.550 


0.521 


6.05(-2) 


2.97 


2.24 














0.528 


0.498 


6.08(-2) 


2.98 


2.24 X 


0.500 


0.521 


6.30( 


-2) 


3.00 


2.22 


0.500 


0.471 


6.09(-2) 


2.98 


2.24 


0.497 


0.518 


6.31( 


-2) 


3.00 


2.22 X 












0.450 


0.469 


6.29( 


-2) 


3.02 


2.22 


0.450 


0.422 


6.02(-2) 


2.99 


2.25 


0.400 


0.416 


6.13( 


-2) 


3.02 


2.24 


0.400 


0.374 


5.82(-2) 


3.00 


2.27 


0.350 


0.362 


5.82( 


-2) 


3.03 


2.27 


0.350 


0.327 


5.50(-2) 


3.00 


2.30 


0.300 


0.309 


5.34( 


-2) 


3.02 


2.32 


0.300 


0.281 


5.04(-2) 


3.00 


2.35 


0.250 


0.257 


4.69( 


-2) 


3.02 


2.39 


0.250 


0.236 


4.44(-2) 


3.00 


2.42 


0.200 


0.204 


3.88( 


-2) 


3.01 


2.47 


0.200 


0.190 


3.69(-2) 


3.00 


2.53 
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p = 0.1 



a/mi = 8 



Roche Sequences 


Irrotational Roche-Riemann Sequences 




02/01 




J 


R/rs 


03/01 


02/01 


02 


J 


R/rs 


0.950 


0.961 


1.23(-2) 


2.82 


4.61 


0.950 


0.949 


1.53(-2) 


2.78 


4.34 


0.900 


0.919 


2.34(-2) 


2.75 


3.88 


0.900 


0.895 


2.80(-2) 


2.73 


3.70 


0.850 


0.874 


3.35(-2) 


2.73 


3.53 


0.850 


0.840 


3.86(-2) 


2.72 


3.41 












0.849 


0.840 


3.87(-2) 


2.72 


3.41 t 


0.832 


0.858 


3.69(-2) 


2.73 


3.45 t 












0.800 


0.828 


4.27(-2) 


2.74 


3.32 


0.800 


0.785 


4.76(-2) 


2.72 


3.24 


0.750 


0.781 


5.09(-2) 


2.74 


3.19 


0.750 


0.729 


5.49(-2) 


2.73 


3.13 


0.700 


0.732 


5.80(-2) 


2.75 


3.09 


0.700 


0.674 


6.07(-2) 


2.74 


3.06 


0.650 


0.681 


6.39(-2) 


2.77 


3.02 


0.650 


0.620 


6.51(-2) 


2.75 


3.01 


0.600 


0.630 


6.85(-2) 


2.78 


2.97 


0.600 


0.568 


6.82(-2) 


2.77 


2.98 


0.550 


0.578 


7.17(-2) 


2.79 


2.95 


0.550 


0.516 


6.99(-2) 


2.78 


2.97 












0.529 


0.496 


7.02(-2) 


2.78 


2.97 X 


0.500 


0.525 


7.33(-2) 


2.81 


2.94 


0.500 


0.466 


7.03(-2) 


2.79 


2.97 


0.494 


0.519 


7.33(-2) 


2.81 


2.94 X 












0.450 


0.472 


7.32(-2) 


2.82 


2.94 


0.450 


0.418 


6.93(-2) 


2.81 


2.98 


0.400 


0.418 


7.13(-2) 


2.84 


2.97 


0.400 


0.371 


6.69(-2) 


2.82 


3.02 


0.350 


0.364 


6.75(-2) 


2.85 


3.02 


0.350 


0.324 


6.31(-2) 


2.84 


3.07 


0.300 


0.311 


6.17(-2) 


2.87 


3.09 


0.300 


0.279 


5.76(-2) 


2.86 


3.14 


0.250 


0.258 


5.39(-2) 


2.90 


3.21 


0.250 


0.234 


5.05(-2) 


2.89 


3.26 


0.200 


0.205 


4.41(-2) 


2.93 


3.38 


0.200 


0.189 


4.17(-2) 


2.93 


3.44 
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P 


0.1 


1 


[7]'s results 


r I SCO tot 


6.5 


7.0 


Our results 


Rfscol^tot 


6.6 


10.4 


Our results 


TJIRRE I Ttyf 

^ISCO 1 ^^^tot 


6.6 


9.9 



TABLE VII. 



p = 1 





a/mi 


5 


8 


[18] 's results 


Tm/Mtot 


8.5 


12 


Our results 


Rfscol^tot 


10.4 


14.5 


Our results 


RlSCo/^tot 


9.9 


14.0 



TABLE VIII. 
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FIGURE CAPTIONS 



Fig. 1(a). The effective potentials of a test particle in the Schwarzschild metric (left) 
and the pseudo-Newtonian potential (right) as a function of the normalized distance R/rg- 

The vertical axes denote ^sch = ^ {l - ^) + l) - 1 and ^pseudo = + 

respectively. The dots express the place of circular orbits. The value of the effective potential 
of the Schwarzschild black hole at the infinity is shifted to zero. 

Fig. 1(b). The fractional deviation of the circular orbits of the pseudo-Newtonian poten- 
tial from those of the Schwarzschild effective potential as a function of the angular momen- 
tum. Circular orbits with different angular momentum related as Jpseudo = {^Y^'^Jsch are 
compared (see text). 

Fig. 2(a). The total angular momentum J = Jtot/{'^i'^2(^s/^toi)^^^} of the equilibrium 
sequence as a function of R/vg in the case that the mass ratio is p = 1. Thin and thick lines 
express the REs and the IRREs cases, respectively. Solid, dotted, and dashed lines denote 
a/rrii = 3, 5 and 8, respectively. 

Fig. 2(b). The axial ratios 02/^1 and a^/ai of the equilibrium sequence as a function of 
R/ts in the case of the mass ratio p = 1. The top figure represents 02/01, and the bottom 
does 03/01. The conventions are the same as in Fig. 2 (a). 

Fig.3(a). The total angular momentum J of the equilibrium sequence as a function of 

R/ts in the case of p = 0.1. The region around the minimum of the lines is magnified and 
shown in the upper right corner. The conventions are the same as in Fig. 2 (a). 

Fig. 3(b). The axial ratios 02/ai and 03/01 of the equilibrium sequence as a function of 
R/ts in the case of p = 0.1. The top figure represents 02/01, and the bottom does 03/01. 
The conventions are the same as in Fig.2(a). 

Fig. 4. The relation between the mean radius of the primary o/mi and the separation 
of the binary Ri scolds at the ISCO in the case of p = 1. Thin hues denote the case of 
the Newtonian potential as an interaction potential, and thick hues denote the case of the 
generalized pseudo-Newtonian potential. The solid and dotted lines represent the IRREs 
and the REs, respectively. 

Fig. 5 The relation between the mean radius of the primary a /mi and the separation of 
the binary Ri scolds at the ISCO in the case of p = 0.1. The conventions are the same as 
in Fig.4. 
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